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easy to handle with CAD/CAM programs without any conversion or approximation. 
The representation was applied to an existing transonic fan geometry to investigate 
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on the aerodynamic performance, and the associated effect on centrifugal stress levels 
was obtained. The investigation demonstrated the suitability of the package to be 
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i INTRODUCTION 


One of the most challenging components to design in an aircraft engine is the 
fan. or the compressor. The compression system has to supply the engine with the 
required mass flow rate, pressure ratio and stall margin with the least possible losses, 
weight and frontal area. 

To satisfy the compressor design requirements, engineers have developed com- 
plex design systems and analysis tools. However, with the conventional design meth- 
ods being pushed to the limit, researchers have started to investigate new concepts 
that were not easy to investigate in the absence of the current computational tools 
several years ago. 

One of the concepts receiving considerable attention currently is sweep. It 
was introduced into transonic fan design to reduce the shock loss. Several early 
sweep attempts were reported, which had little success. More recent attempts have 
shown promise. However, the problem is not fully understood, and a clear isolation 
of sweep effects is not yet available. 

Viewing sweep as a promising concept, an attempt was made here to investi- 
gate the incorporation of sweep into the design process. The current design system 
was reviewed first to determine the best way to approach the problem. From that 
review it was concluded that the current design system is not suitable for the problem 
at hand specifically in the way that the geometry is handled. The geometry used in 
current design systems is a 2D representation of blade sections that can be moved cir- 
cumferentially or axially in the structural design phase. This would be unacceptable 
when designing with sweep. Moving the blade sections circumferentially after decid- 
ing on a sweep schedule will alter the flow significantly. The transonic flow through 
a blade row with sweep is essentially three dimensional. This led to the conclusion 
that a 3D blade geometry representation is badly needed if a rational 3D blade design 


procedure is to be developed successfully. 


Bezier surface representation was found to be suitable for blade representation. 
Having few control parameters and offering flexible manipulation, it was viewed as a 
good choice for 3D blade geometry representation, one that might lead to 3D blade 
design optimization. The other significant advantage of the Bezier representation, 
besides being compact, is that it is read easily by CAD/CAM programs with no 
approximation. 

The representation was packaged in a form suitable for turbomachinery blad- 
ing and compatible with the existing design system. A test case was applied to prove 
the concept and determine the suitability of the package to the design problem. The 
representation was found to be very flexible and capable of rational geometry manip- 
ulation. It was also found that the blade design parameters used currently can be 
represented in terms of the Bezier parameters. 

The effect of sweep on blade row performance was investigated. The sweep 
was implemented in a manner that kept the blade parameters unchanged (maximum 
thickness, maximum thickness location, section profile, etc.). This was in contrast to 
previously reported sweep attempts that had to accept other changes to the design to 
meet steady and unsteady stress constraints. The geometry represented by the Bezier 
format was imported into IDEAS with no approximation, and a purely centrifugal 
stress analysis was conducted. 

The results of the test cases were that the swept-back rotor efficiency was 
better than the swept-forward, and the rotor pumped higher mass flow rate at the 
same speed and pressure ratio. The effects of sweep were qualitatively similar when 
changes were made to the leading edge within the grid generation software.. 

The geometry package still needs further development and the tools supporting 
it also need to be developed to be compatible (grid generation in particular). The 
obtained results are viewed as a proof of concept and not as a definitive study of the 


effects of sweep. 


mle PROBLEM REVIEW 


Aircraft compressor design is a very complicated process; one that involves 
aerodynamic, structural, manufacturing and maintenance disciplines. Each of these 
introduces very tight constraints that limit the design space. In addition, the tools 
that are used in each discipline dictate what information is needed from the others. 
The ease of interfacing and data sharing is very critical. Usually, to go from one engi- 
neering or manufacturing group to another, a conversion and approximation process 
takes place. 

The aerodynamic design is necessarily limited to what analytical tools are 
available. Typically, the designer arrives at a geometry using analytical methods 
based on simplifed models. He then verifies that geometry by analysis using avail- 
able computational tools. Design methods over the past decades have evolved from 
the use of very simple one-dimensional mean line models, to the incorporation of 
"simple radial equilibrium”, to quasi-three dimensional (Q3D) methods for arbitrary 
flow paths. The quasi three-dimensional methods became practical when methods 
such as stream line-curvature for calculating the axi-symmetric through flow became 
available. 

Similarly, the computational analysis tools have evolved from two-dimensional 
potential plus integral boundary layer, to three-dimensional Euler plus differential 
boundary layer, to fully three-dimensional Navier-Stokes methods. The proven capa- 
bility of computational fluid dynamics (CFD) to reproduce complex three dimensional 
flow fields has recencly led to the use of expressions such as "CFD design methods” 
and, even more recently, to the idea of trying to rationally optimize blade design to 
include ”sweep” and ”lean” as well as profile shape. In order to put the present work 
in perspective, a review is given of the classical Q3D and recent CFD design methods, 


and of published attempts to introduce sweep. 


A. DESIGN METHODS 
1. Quasi Three Dimensional Method 


In 1952, Wu [Ref. 1] developed a general theory of three-dimensional flow 
in subsonic and supersonic turbomachines. It was for steady nonviscous rotational 
and irrotational flows. There was no limitation on whether the machine should be 
axial, radial or even of mixed flow type. The solutions of the direct and inverse 
problems were possible through the definition of a specific combination of relative 
stream surfaces. Referring to Wu’s illustration shown in Fig. 1, solving for the two- 
dimensional flow on two stream surfaces identified as $1 and $2 equivalently solves the 
three-dimensional problem. S1 represents the ’tangential-axial’ stream surface while 
S2 represents the ’radial-axial’ plane. In general, these surfaces can be of arbitrary 
shape to satisfy the equation of motion. The flow channel can have as many 51 and 
S2 surfaces as the problem may require. An iterative procedure is used between the 
two surfaces to solve for the flow field. The solution on each surface updates the flow 
on the other surface. 

Since it was introduced by Wu, the 51 and $2 method has been applied in 
different ways. Simplifications were introduced which varied from one method to 
another. Most of the methods used one S2 surface and several 51 surfaces. The most 
widely used simplification is the assumption that the $1 plane is axisymmetric. The 
assumption of linearly varying streamtube thickness and radial position allowed the 
omission of in-blade computing stations. On the other hand, an extension to the 
theory to handle rotational transonic flow has been reported [Ref. 2]. The use of a 
number of $1 and S2 surfaces was also reported [Ref. 3]. 

The quasi-3D design approach [Ref. 4, 5, 6] is based on the $1 and S2 formula- 
tion, using the assumption of axisymmetric $1 surfaces and axisymmetric conditions 
between blade rows. First, given a mass flow rate, rotational speed and work, the 
annulus height and mean-line flow angles are determined using a one-dimensional 


mean-line method. Then a through-flow analysis is used to determine the span-wise 





Figure 1. Intersecting $,52 Surfaces in a Blade Row [Ref. 1] 


variation of inlet and exit flow angles to the blading and the streamtube thickness 
variation. The blade section geometry, which is defined on the axisymmetric stream 
surfaces, is produced using blade-to-blade inverse or direct methods, or using empir- 
ical correlations. The design system is not fully analytic; it has to be complimented 
with experience and empiricism. The selection of loss coefficients, incidence angles 
and diffusion factors is normally based on correlations derived from experimental 
results. Such data are updated when test results for newly designed blading are 
available. Most design data are proprietary and not available in the open literature. 


In 1985, Jennions [Ref. 6] presented a Q3D blade design method that relied 


on a Rolls-Royce database. The system consists of a through-flow program, a blade- 
to-blade program and a stacking-line program. Illustrated in Fig 2, the through-flow 
program is based on the streamline curvature method. The flow channel is divided 
into computing stations from inlet to outlet including the blade region. The radial 
equilibrium equation is solved on each station to refine the streamline path. The radial 
pressure gradient calculated from radial equilibrium is affected by the streamline 
curvature in the meridional plane, centrifuging of the fluid, axial pressure gradient, 
blade force, dissipative body force and perturbation (non-axisymmetric) effects. The 
flow is treated as axisymmetric in the duct areas while in the blade area a passage 
averaging is used to avoid any indeterminacy in defining the $2 surface. At each 
computing station, initial mid-span static pressure is prescribed to calculate the mass 
flow rate through the station. Accordingly, the static pressure is adjusted to achieve 
the desired mass flow. 

The blade-to-blade program uses the stream surface geometry from the through- 
flow analysis to determine the blade section shape. The design is made either in a 
direct or inverse mode. The program satisfies the aero and structural constraints for 
the designed blade section. 

The stacking line program stacks the designed blade sections from hub to tip 
to form the three-d:mensional blade geometry. It stacks the sections at any circum- 
ferential position in the annulus as long as it is structurally acceptable. The reason 
for this freedom 1s that the blade sections are viewed as isolated from each other. 

While relatively easy and fast, the Q3D method has significant limitations. 
The particular approach described by Jennions [Ref. 6] is inviscid and does not 
account for blade lean or sweep. Also secondary flows are not predicted and the flow 


on each blade section does not communicate with the rest of the flow. 


Z. Three Dimensional Method 


Fully three dimensional CFD methods have been incorporated into the design 
process recently due to the limitations in the Q3D approach. Whether the method 
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Figure 2. Streamline Curvature Computing Stations and Stream Surface [Ref. 5] 


is based on Navier-Stokes equations or Euler equations with boundary layer correc- 
tions, the new approach starts with a geometry arrived at by stacking blade sections 
designed with Q3D methods. The sections may or may not have been ’optimized’ 
for the assumed axisymmetric flow. Then the designer applies the three dimensional 
solver to the geometry and examines the flow field solution. The blade sections are 
then changed and restacked to give a new 3D geometry. The process is repeated until 
the desired performance, both on and off-design, is achieved. 


Examples of using 3D CFD codes in the design process for transonic compres- 


sor stages are seen in the work of Denton [Ref. 4] and Sanger [Ref. 7]. In both cases 
a 3D Euler code, incorporating a viscous boundary layer by modelling, was used to 
analyze the flow through a geometry that was the result of stacked blade-element 
sections. In the process of modifying the geometry, changes were made to parameters 
defining individual elements, and the code was then reapplied to analyze the new 
blade. The geometry packages were similar in that blade elements were defined by a 
camber line, a distribution of thickness, and leading and trailing edge shapes. The 


elements were then stacked into a blade by defining the stacking line. 


B. SWEEP ATTEMPTS 


In external aerodynamics, swept wings are used very successfully to reduce the 
drag associated with shock waves. Extending the sweep concept to turbomachinery 
is not a simple extension however. Conditions are not uniform along the blade span 
and the internal shock structure is completely different from that found in external 
aerodynamics. 

While the relations between free-stream velocity, wing sections and shock wave 
structure is easily determined in external aerodynamics; the opposite is true in tur- 
bomachinery. Each blade section operates at a different relative flow velocity and 
flow angle. Centripetal force, tip leakage, secondary flows and non-uniform inlet flow 
conditions contribute to the problem’s complexity. What is commonly observed in 
transonic rotors is a contained "passage shock”. Based on this observation a loss 
model was developed by Miller et al [Ref. 8], which was later modified by Wenner- 
strom [Ref. 9]. 

The initial model shown in Fig. 3, assumed a normal shock in the passage 
extending from point A to point B. The blade spacing affects the location of point 
B which is determined by a perpendicular to the mean passage camber line going 
through the following blade leading edge point A. The shock loss is found by averaging 


the free-stream Mach number and the Mach number resulting from expanding the flow 


around the leading edge to point B. A more detailed model was presented in the same 
paper (see Fig. 4) that accounts for the varying losses across the passage due to the 
varying Mach number. The variation in Mach number is calculated by accounting 
for the bow shock and the expansion waves across the channel. The compressor loss 


1s computed by adding the losses for small stream tubes along the blade span. 
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Figure 3. Passage Shock-Wave Configuration for Estimating Shock-Loss Levels [Ref. 
8] 


BOW WAVE 


SUPERSONIC 
— 

UPSTREAM 

FLOW 





~-——tt some Une 


Figure 4. Flow at Inlet of Cascade of Airfoils Operated with Supersonic Relative Inlet 
Mach Number [Ref. 8] 


Calculating the losses with the described model resulted in shock losses higher 
than the total losses obtained from experimental results. The consistent loss over- 
prediction was due to the two-dimensional treatment of the shock [Ref. 9]. The 
second model treated that limitation by considering the shock obliquity in the third 


dimension Fig. 5. The model follows the first model in the calculation of stream Mach 


number and determining the shock location on the suction surface. The shock surface 
is then constructed and the relative angle between the surface and the free stream 1s 


calculated on each section. 


SHOCK SURFACE 





Figure 5. Three Dimensional Shock Surface Between Two Blades [Ref. 9] 


1. NASA QF-12 


In an effort to reduce acoustic noise through the reduction of shock wave- 
generated pure tones, NASA sponsored the design of a swept fan stage [Ref. 10] with 
relative inlet Mach number of 1.588. The compound sweep reduced the relative Mach 
number to an effective Mach number of 0.91 at the mean and tip with a hub Mach 
number of .83. The fan did not meet the aerodynamic design goals and the acoustic 


measurements, although below expectations, showed some improvements. 


2. Neubert et al. 

In 1986, two transonic fans were designed [Ref. 10], one was swept back 
and the other was conventional for the purpose of comparison. The swept design was 
designed to be shock free with a pressure ratio of 2.21, tip speed of 1558 ft/s corrected 
flow of 1.51b/s, relative Mach number of 1.6 and adiabatic efficiency of 88.1%. An 


estimated gain in efficiency of 1.5% was estimated. The loss reduction was due to the 
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elimination of the calculated shock wave loss and the loss due to the shock boundary 
layer interaction. 

The rotor was designed using a 3D Euler solver with a two-dimensional integral 
boundary layer procedure to estimate the viscous effects. The design procedure is 
shown in Figure 6. The numerical solution showed that the Mach number contours 
on the swept suction surface was spaced out, suggesting that the shock disappeared 
over most of the blade span. The tip section had some concentration of Mach number 


contours, but less than the conventional rotor. 
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Figure 6. Swept Fan Design Procedure. [Ref. 10] 


The test results for the two rotors showed that the swept rotor had poorer 
performance than the unswept rotor. The pressure ratio was 0.2 lower, the cor- 
rected mass flow was less by 0.9/b/s and the stall margin was reduced substantially 
to (0 — 6)% compared to 13%. There was an improvement of 0.1% in corrected 
polytropic efficiency (corrected to match relative Mach number). The difference in 


performance between the swept and unswept rotors was related to the differences 


1a) 


found in span-wise loading and mass fraction distribution. Detailed measurements 


showed the existence of a shock surface in the outer half of the swept rotor. 


3.  Puterbaugh et al. 

Puterbaugh et al [Ref. 11] designed a swept back rotor derived from a base 
line geometry [Ref. 12}. The swept design was mainly to evaluate the use of the 
3D-shock loss model in the design system, and the overall effect of blade sweep on 
compressor performance. The detailed design (started from the base line geometry) 
involved optimizing the blade sections to obtain a circumferentially-averaged static 
pressure distribution that rose linearly from the leading edge to about 3/4 chord 
and decreased smoothly to give zero slope at the trailing edge. Relative flow angles, 
incidence and deviation angles were used to derive the general airfoil shapes. The 
design parameter (incidence, deviation, loss, work distribution, blockage) selection 
was based on experience and previously-published designs. The total pressure loss 
was obtained from the diffusion factor correlation combined with the 3D-shock loss 
model. Then sweep was added to reduce the shock loss in the rotor. The amount of 
sweep was increased radially as the Mach number increased to obtain a leading edge 
profile optimized for minimum shock loss. The structural design was considered while 
deciding on sweep schedule to ensure structural integrity. 

The rotor design goal was for specific flow rate 214.7kg/sec/m’, pressure ratio 
2.04, efficiency 94.0%, tip speed 457.2m/sec. The test results showed that the flow 
rate was achieved while the pressure ratio fell short by 0.07, efficiency was lower by 
3% and the stall margin was 11%. The shock model was found to underestimate 
the shock sweep angle at mid span while overestimating the angle at the tip. This 


resulted in lower mid-span work addition. 


an Wadia et al. 


In an effort to better understand sweep, Wadia et al [Ref. 13] carried out an in- 


vestigation that included three swept rotors. They designed two swept backward and 
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one swept forward. The swept-forward rotor and one of the swept-back rotors were 
swept through tangential lean and required barreling the pitch chord (pitch chord is 
the midspan chord [Ref. 14]). The other swept-back rotor was swept through tan- 
gential lean only. The designs had the same tip speed, inlet-corrected flow, pressure 
ratio and passage area distribution. The chords and solidity were changed due to me- 
chanical constraints. In addition, the maximum thickness and location of maximum 
thickness could not be held fixed and were viewed as a part of the sweep package. The 
forward sweep required added root thickness due to increased blade length, which af- 
fected the blade flexural frequency. The two barreled rotor weights were significantly 
increased. The forward-swept rotor required a change of material to one of superior 
mean stress capability. 

The rotor experimental results were compared to each other and to the base- 
line rotor (the same as in Puterbaugh et al). The barreled, swept-back rotor showed 
the worst stall margin (of only 4% compared with 9% for the baseline) with a nearly 
unchanged (even slightly increased) peak efficiency. The unbarreled swept-back rotor 
stall margin was 7% and showed an increase in the peak efficiency of 0.9%. The 
forward-swept rotor had the highest stall margin at 13.5% with an unchanged ef- 
ficiency. The experiment was supported with CFD analysis which showed that the 
boundary layer accumulated in the tip region for the swept-back rotors. This accumu- 
lation was identified as the main cause of stall margin reduction. The forward-swept 
rotor was found to have relatively high mass flow rate all the way down to lower 
speeds. This was interpreted to be due to the unique characteristic of the swept- 


forward blade sections to untwist (open) as as the rotational speed ws reduced. 


C. LIMITATION IN CURRENT DESIGN SYSTEM 


There has been a major effort made in recent years by the turbomachinery 
community to improve compressor design. This is evidenced by the increasing number 


of designs derived using 3D solvers, including a number of swept rotors. Although the 
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effort is significant there have been no reports of designs using Navier-Stokes solvers. 
This may be attributed to the fact that they were not available until recently and 
require large execution times. 

One of the main limitations in current design systems is the method for de- 
scribing the blade geometry. The current methods were quite satisfactory for con- 
ventional Q3D designs, before the advent of capable 3D analytical tools. With the 
introduction of sweep and general blade shapes, designing the blade as independent 
two-dimensional sections and then stacking them can lead to several problems. Due 
to the absence of constraints in stacking the blade sections there can be waviness in 
the final blade surface which is costly to machine. In addition, satisfying geometrical 
conditions of continuity (of contour, slope and curvature) in each two-dimensional ge- 
ometry does not necessarily guarantee that condition in the third dimension. In fact, 
current aerodynamic design blade geometry packages communicate neither well with 
the structural analysis programs nor with the CNC manufacturing machines. The 
aerodynamic shape is only approximated in the manufacturing process. Finally, with 
current geometry packages, the full 3D blade shape optimization problem is unthink- 
able. For a blade to have from 11 to 22 sections along the span and around twenty 
coefficients for each section, the optimizer would have from 220 to 440 parameters to 
perturb. In view of the time it takes the solver to complete one cycle, it is obvious 


that there is no hope for an optimized solution. 


D. WHAT NEEDS TO BE DONE 


The geometry packages currently used to design turbomachinery blading are 
based on the stacking of 2D cascade blade elements. While this representation is 
required in the Q3D design process, it is not required in the 3D analysis process. 
Moreover the continued reliance on this representation is inhibiting the generation of 
designs which are truly optimized on the basis of aero, structural and manufacturing 


considerations. What is therefore required is a geometry representation that 
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1. is fully three dimensional in character, 
2. involves fewer parameters than current representations, 


3. 1s compatible with structural analysis and CNC machine control programs, 
without additional approximations, 


4. can represent existing geometries to acceptable accuracy, 


5. can be programmed to allow changes to the geometry to be made easily by 
the designer, and those changes to be related to Q3D blade parameters. 


The development of such a package and its application to the introduction of sweep 


into a rotor design, was the goal of the present work. 
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III. GEOMETRY REPRESENTATION 


A. CHOICE OF GEOMETRY REPRESENTATION 

Objects can be defined in many different ways. A general curved shape may 
be represented using a great number of small straight lines or with combinations of 
several arcs. Many years ago, parts were shaped by deforming a flexible strip of 
metal by hanging weights at specified points to form a smooth curve. That flexible 
strip was called a spline. Today, with major advances in manufacturing and in the 
accuracy that can be achieved in the machining of a part described with numerical 


data, parametric representation has been adopted for representing spline shapes. 


1. Different Geometrical Representations 

Both parametric and functional (non-parametric) representation can be used 
to describe curves and surfaces. Usually the parametric representation is more con- 
venient and more flexible. A circle for example can be represented in non-parametric 


form by the single equation 
(x —a)*+(y—b) =r’ (III.1) 
or in a parametric form, with a variable parameter 0, by the two equations 


s = ?Tcos(O)mme (IIT.2a) 
rsin(@) + ( EZ 15) 


ce 
| 


The parametric representation can be of different forms. Some of the repre- 
sentations available are Hermite, Bezier, and uniform B-splines which are described 


respectively by 


Q(t) = (2t° — 3t? +1)P, + (—22° 4 327) P, + 
(t? — 20? + t)R, + (t —t’)Rg (IIT.3) 


ef 


Q(t) = (1—2)P; + 3t(1 —t)P, + 307(1 — 2) Pa + PY (III.4) 





1—t) 3t° — 6t? + 4 
gi) = | ; Pig + et 
__ 243 2 t? 


Q is a coordinate of a point on the curve. The P’s and R’s are the corresponding 
blending function weights. 

Hermite curves are defined by two end control points and the tangents at the 
end points. Fig. 7 shows the blending functions used to generate Hermite curves (the 
polynomials which are coefficients of P,, Ps, Ry, and R, in Eq. III.3) and a Hermite 


CUTIVE. 
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Figure 7. Hermite Representation 


Bezier curves are defined by four control points, two of which are the end 
points. The other two are not on the curve but with the end points define the 
tangents at the end points. Fig. 8 shows the Bernstein polynomials (which are the 
polynomial coefficients of the P’s in Eq. III.4) and a Bezier curve with its four control 


points. 
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Figure 8. Bezier Curve Representation 


The B-spline is a general representation, which does not intersect with any of 
its control points. It defines multiple segments of curves with a uniform parameter. 
Fig. 9 shows the basis functions and a one segment B-spline. 
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Figure 9. B-Spline Representation 
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2. Advantages and Disadvantages 

Each of the above-mentioned representations has particular advantages and 
disadvantages. Parametric representation is advantageous to the functional represen- 
tation in the ability to represent shapes easily, define a part of a standard geometry, 
and to represent a complex geometry. All have defined tangents. More complex 
geometries may be represented by higher-order functional representations which are 
hard to control between points, and can wiggle. The polynomial coefficients have no 
direct link to the geometry shape. The curve or surface functional representation has 
to be complemented with the range of definition. The advantage of the functional 
representation is the ease of determining if a given point lies on the curve or not. 

Between the different parametric representations, there are advantages and 
disadvantages although there is a transformation matrix between the corresponding 
forms. Table I shows a comparison between the properties of the different forms. 


Items appearing in the table will be addressed in the following discussion. 


3. Selected Representation 

The representation selected here for blade geometry is the cubic Bezier repre- 
sentation. In addition to the properties that will be presented in the following sections, 
they hold the particular advantage that the control points are the end points (lead- 
ing and trailing edge points) and points that control the slopes at those end points. 
They can be imported into any CAD/CAM program with no approximation as they 
are a special case of NURBS (Non Uniform B-Splines). They are widely used for 


representing curves and surfaces. 


B. BEZIER CURVES 
1. Mathematical Representation 


Bezier curves can be represented in many different forms. The forms are 
equivalent and give the same results. Each variable is defined by its own control 


points and the same blending functions. For curves of dimensions higher than one, 


20 


Ss a _ Bezier | Uniform | Nonuniform 
Pe | Spline | Spine 

Convex hull defined 

by control points 


Interpolates some 
control points 


Interpolates all 
control points 


Ease of 
subdivision 


Continuities inherent 
in representation 


Continuities 
easily achieved 


Number of parameters 
controlling a curve 
segment 





Table I. Comparison of Properties of Different Forms of Parametric Cubic Curves 
[Ref. 15]. 


the control points represent vectors in the corresponding space. All the following 


representations will be given for the cubic form. 
ae Polynomial Form 


This form of the Bezier representation is not commonly written but is 
included here because polynomial representations have been the most widely used in 


generating blade shapes: 
Q(t) = (—po + 3p: — 3p2 + ps)t* + (3p — 6p: + 3p2)t” + (—3po + 3p1)t + po (IIL.6) 


In the above equation, the constant term, po, is the first control point and the last 


control point, p3, is sum of the polynomial coefficients. 


ah 


b. Bernstein Form 
The Bernstein polynomials are recursive polynomials. The higher order 


polynomial is calculated from the next-lower degree polynomial, namely 


Bee (1 — ee (7 eee ee (III.7) 


The Bernstein polynomial is given in the binomial form as 


it, Ne | 
aE )iele Salant). (II1.8) 
1 
The Bezier formula is then . 
Q(t) = > piB?(t) (III.9) 
2=0 
Ce Matrix Form 
The matrix form of the Bezier representation 1s 
1-3 3 +l te 
0 3 -6 3 Ee 
Q(t) = [po Pi P2 ps! (IIT.10) 
00 3 -3 i 
pO 0 80 Lap rE | 
The matrix form is useful for computer implementation, especially if matrix multi- 


plication is hard wired. 


2. Cross Plots 
The Bezier curve in the zy plane or the ryz space is composed by cross plotting 
the coordinates in terms of each other for given values of the parameter t. This 1s 


illustrated in Fig. 10 for the two-dimensional case. 


3. Effect of Moving Control Points 

Bezier curves are sensitive to control point movements. For example, moving 
the third control point of the Bezier curve shown in Fig. 11 by 20% of the tangent 
length in the positive and negative directions changes the shape of the solid curve to 
become the broken curves. Notice the curve is moved in the same direction as the 


control point. 
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Figure 10. Cross Plot 


4. Properties of Bezier Curves 

Bezier curves have properties that make them particularly attractive for design 
and graphics. ‘They are invariant under translation, rotation and scaling but not under 
perspective projection. Convex hull’, linear combination and end point interpolation 


are also interesting properties. 
a. Convex Hull 


The ’convex hull’ property of the Bezier curve is that the curve is 
contained inside the polygon of the control points defining the curve, as shown in 


Fig. 12. This is helpful in knowing the extent of the curve without computing it. 
5. Invariance 


The Bernstein polynomials are afhne map invariant, which is stated as 
Wr-o BY (t) = 1. This makes the Bezier curves invariant under translation, rotation 
and scaling. That allows applying transformations to the four control points and then 
calculating the transformed curve. The equation 
ON > pi Br (t) +a= 2 (p; + a) Br (t) eins 
j3=0 


j=0 


shows that translation applied to the control points, or to the curve, 1s equivalent. 
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Figure 11. Control Point Movement Along The Tangent Direction 
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Figure 12. Convex Hull Property 


The curves are also invariant under linear combination, since 


S| (ab; + Bc;)Ba(t) =a >> b,BP(t) + BY — c; Br(t) (IIT.12) 
j=0 j=0 j~0 
GC: Location of Mazimum Effect 


Examining the Bernstein polynomials, each one is seen to have only 


one maximum. That maximum is at t = + [Ref. 16], where z is the number of the 


a 
n 


term in eqn. III.4 and n is the order. So to change the curve at one third of the 


parameter, the second control point movement would give the maximum effect. 


24 


5. Differentiation 

The derivative for a Bezier curve can be obtained simply by differentiating any 
of the given forms. It is worth mentioning though that the derivatives have direct 
relationships to the control points. By looking at the polynomial form in Eq. III.7 
it can be seen that the order of the polynomial for the derivative will be reduced by 
one and that the derivative at zero will be three times the difference between the first 
two control points. This is not only true at the end points. The polynomial for the 
derivative for the whole curve is a Bezier curve of order one degree less than the curve 
itself (second order in our example). The control points for that curve are simply the 


difference between the successive central points of the curve being differentiated, or 


n-1 
HO) = 9 Oi — 2) BIO (111.13) 

6. Continuity Between Curves 

CAD/CAM application requires more than one curve to be joined together. 
The requirement on connecting two curves may vary from one case to another. An 
application may require curves joining with first derivative being continuous or be 
both first and second-derivative continuous. For first-derivative continuity, the two 
curves have to have equal tangents, or one has to be a multiple of the other. That 
would be satisfied by forcing the three joining control points to be collinear as shown 
Terie. 13. 

Curvature continuity can be ensured by arranging the control points on each 


side to maintain a common ratio between the segment lengths shown in Fig. 14 


We Difference Between C&G 


A function with a continuous second derivative would be called C? continuous. 
This notion is still valid for parametric representation, with reference to variation in 
the parameter. When cross plotting two parametric curves however, the resulting 
curve might not have the same degree of continuity that was inherent in the original 


curves. In this case we have another type of continuity which is called geometric 
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3 Common Tangent 







bi 


Bezier Cu 


Bezier Curve 1 


i an 


. 
Figure 13. Two Bezier Curves with First-Derivative Continuity 


continuity (e.g. G* for geometric curvature continuity). Having one type of continuity 
does not guarantee the other type of continuity. The geometric continuity is less 


constraining and the more important to design applications. 


8. Subdivision 
One of the interesting properties of Bezier curves is subdivision. Given a Bezier 


curve with a set of control points, we can always find two sets or more of control points. 





ai/b1 = a2/b2 = a3/b3 


Figure 14. Condition for Second-Derivative Continuity for Bezier Curves 
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Fach of these would represent the corresponding segment exactly while maintaining 
C* continuity. The procedure is relatively easy and straightforward. Consider the 
control polygon defined by the control points starting from the first point to the last 
(P0, P1, P2, P3) as shown in Fig. 15 and a parameter t, where subdivision is to occur. 


The control points (a,d, f,e,c) can be found using the equations 






Point of subdivision 


i P3 
Figure 15. Bezier Curve Subdivision 
@ = 7,P1 4 (1 —7.)P0 (III.14) 
b = t,P2+(1-—t,)Pl (III.15) 
ec = t,P3+(1-t,)P2 (III.16) 
d = t,a+(1—-t,)b Cet 7) 
e = t,b+(1—-t;)e (III.18) 
f = t.d+(1-t,)e (III.19) 


We obtain two Bezier curves defined by (P0,a,d, f) and (f,e,¢, P3). 


Zt 


C. BEZIER SURFACES 

Bezier surfaces are defined by sixteen control points. The edges by themselves 
are Bezier curves. The surfaces possess the same properties as the curves. An individ- 
ual Bezier surface is called surface patch. Surfaces can be formed by joining surface 
patches together with the required degree of continuity. Usually surface continuity is 


more difficult to satisfy than continuity for curves. 


1. Mathematical Representation 
The Bezier surface is defined with two parametric variables u,v. The form of 


the Bezier patch is 


Q(uv) = OV py NP (uw) MM(v) (111.20) 


1=0:'7=0 
where V,M .... are Bernstein polynomials. A Bezier surface with its control net is 


shown in Fig. 16. 





Figure 16. Bezier Surface with its Control Net 
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2. Effect of Moving Control Points 
The effect of moving one of the inner control points for the surface shown in 
Fig. 16 is illustrated in Fig. 17. Only one coordinate of the control point was changed 


in this example. 





Figure 17. Effect of Surface Control Point Movement 


3. Continuity Between Surfaces 
The conditions for surface continuity are the same as for the curves with the 


extra condition that the four tangents along the connecting edge shown in Fig. 18 


must be such that al/bl = a2/b2 = a3/b3 = a4/b4. 


D. PARAMETERIZATION 


In the previous section we dealt with the problem of calculating a surface 


from a given set of control points. This is not always the case, if we already have the 
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Figure 18. Tangency Condition for Surfaces 


geometry as a set of data points we would like to obtain the set of control points that 
best represents the surface data points. 

The problem is solved by finding the best parametrization for the given curve 
or surface. Hoscheck [Ref. 17] solved the problem by splitting it into two least- 
squares problems. The first is the least-squares problem of finding the control points 
for the surface that best fits the given geometry data points assuming a given param- 
eterization. The second is a nonlinear least-squares problem that solves for the best 
parameterization given the data points and a set of control points. 

The solution procedure described by Hoschek [Ref. 17] is iterative and proceeds 


as follows: 


1. Find an initial parameterization using a suitable method (arc length for ex- 
ample). 


2. Apply least squares to find the approximating control points. 


3. Calculate the error vector D; at each point (see Fig. 19). 
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. Find the parameter correction at each point by calculating the dot product 
of the error vector and the tangent vector at the corresponding point. Then 
scale with the ratio of the parameter length to the total arc length. 


. Repeat until the solution converges. 







© Fitted Points 


* Data Points D{i+3) 


Figure 19. Error Vector D; for Parametrized Data 
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IV. INITIAL GEOMETRY PACKAGE 


Aa a first step. the suitability of the Bezier representation for describing typical 
2D compressor blade section profiles was examined. It was found that Bezier curves 
were suitable and the blade shape was well represented. Second, a complete blade was 
treated as two separate surfaces, pressure and suction surfaces, and the representation 
was again found to be suitable. After representing the pressure and suction surfaces 
with Bezier patches, leading and trailing (Bezier) surfaces were added. This resulted 
in a complete blade represented in a form that was compact (32 control points and 2 
parameter values) and also easy to handle, and easy to import into different programs 
without approximation. Most importantly, the coordinates at any position on the 
blade surfaces could be retrieved with no loss in accuracy (other than that due to 
machine precision), because the shape was analytically described. 

The package was implemented with current compressor design systems in 
mind. The initial package can read an existing geometry specified in meridl2, meridl3 
and general format [Ref. 18]. The package can also read the control points that repre- 
sent the six surfaces or only the ones that represent the pressure and suction surfaces, 
and then add leading and trailing edge surfaces. Also, to facilitate further develop- 
ment and the addition of a graphical user interface, the package was implemented 
using object-oriented programming. 

The blade surfaces are numbered from 1 to 6 as shown in Fig. 20. The control 
points for each surface are numbered from 0 to 15 as shown in Fig. 21. Where surfaces 
are joined, the surface edge with the lower numbered control points is always on the 
left (when moving clockwise around the blade). For example, surface No. 2 (pressure 
side) control points (3,7, 11,15) will match (are the same as) control points (0, 4, 8, 12) 
of surface No. 3. 

The blade geometry is manipulated by moving control points individually or 


in-groups defined in the program. The groups are defined such that one row of control 
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at 


Figure 20. Blade Surface Numbering 


points for the six surfaces can be moved simultaneously. The movements allowed give 
the freedom to shape the geometry with great flexibility. 

The use of the present representation is new to axial turbomachinery and 
by no means is the present work considered to be complete as an exploration of the 
potential of the representation. The present work only demonstrates proof-of-concept. 
The intent was to find a rational approach to design with sweep. In that context, 
basic blade geometry manipulation routines and prescribed parameter variations were 
implemented. Specifically, two blade sweep schedules, blade lean, and individual 
control point manipulations were programmed to be computed relative to a reference 


geometry. 


A. APPROXIMATING BLADE GEOMETRY 

The blade geometry is read in as r,z,6,02 as defined by either meridl3 or 
meridl2. Meridl3 reads in the blade section with 6, being the mean camber line 
tangential coordinate and @) the tangential thickness angle. Meridl2 reads in the 


pressure surface and suction surface tangential coordinates as 6, and 62 respectively. 
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Figure 21. Control Point Numbering 


In either case, the coordinates are transformed to zx, y, z coordinates and the surfaces 
are treated as two ace surfaces by the program. The surfaces are assumed to 
start and end at their defined leading and trailing edge tangency points respectively. 
Each surface is fitted by a Bezier surface separately using the method presented by 
Hoschek [Ref. 17]. The fitting routine starts with initial parameterization that uses 
the cubic root of the arc length. The linear least-squares fitting uses QR factorization, 
while the parameter correction uses the method given by Hoschek [Ref. 17]. 
The equation for the initial parameterization is 


where s is the distance along the curve. The parameter interval is then normalized 
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Woose 
-0.0572 
-0.0044 
-0.0482 
-0.0891 





-0.0777 

-0.0145 

0.0002 

-0.0340 

-0.0673 
0 


Table II. Different Parameterization (t = parameter, e = error; subscripts a,s,c are 
arclength, square root and cubic root). 


by calculating 
dt; 


ed 


Table II shows the accuracy of predicting the initial parameter this way versus using 





Gt (IV.2) 


the arc length or the square root of the arc length. The first three columns of the table 
were generated by calculating data points from Bezier control points with an unevenly 
spaced parameter t. The next six columns compare the various approximations. 

Since Bezier surface corner points are control points, the values for the corner 
points are substituted into the equations to fix the corners. The number of equations 
is then reduced by four. The problem becomes a classical least-squares problem in 12 
unknowns. 

The nonlinear part is solved by calculating the error vector at each point and 
calculating its inner product with the surface tangent with respect to both u, and 
v parameters [Eq.III.20]. The result is scaled with respect to the corresponding arc 
length and added to the parameter value at each corresponding point. 

Two-dimensional blade sections were tested to determine the suitability of the 


representation. ‘T'wo test cases are shown in the Fig 22 and Fig. 23. 
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Figure 22. Double Circular Arc Blade Section [Bezier Approximation] 


B. ADDING LEADING AND TRAILING EDGE SUR- 
FACES 


The two Bezier surfaces that represent the blade pressure surface and suction 
surface need to be closed to form the blade geometry. In defining the leading and 
trailing surfaces there are two considerations. First, the connecting edges have to 


maintain the appropriate continuity requirement. Second, a certain amount of flexi- 





bility should be provided so that the leading edge bluntness can be changed. Joining 
the two surfaces by one surface, with even the least amount of continuity, would not 
satisfy these requirements. Two surfaces were judged able to satisfy curvature con- 
tinuity at the pressure and suction sides while having at least tangent continuity at 
the edge connecting the two added surfaces together. 

The leading surfaces are defined by the constraints on the ratio of the tan- 
gents to the blade surfaces, and their common edge. This determines the chord-wise 
extension added by the leading (and trailing) surfaces. Since the requirement is to 
have constant ratio of tangents along the surface edge [Fig. 17], and the arc length 
is larger at the tip than at the hub, the tip extension will be larger than the hub 
extension. The chord length is therefore adjusted by applying subdivision to bring 
the chord length back to the desired length. 

After adding the leading and trailing surfaces the blade is fully defined, with 
the six surfaces numbered as shown in Fig 21 from 1 to 6 (clockwise starting from 


the edge connecting the two trailing edge surfaces). 


C. EXAMPLES OF GEOMETRY MANIPULATION 


The manipulations implemented in the present package are of two types. One 
is a general point manipulation in which each individual control point can be moved in 
any of the three coordinate directions (r, z,@). Second is a programmed manipulation 
in which a defined set of control points is moved in any of the three directions (r, z, 9). 
This second type includes a manipulation in which the blade is swept along the blade 
chord. This is applied either to the top row of control points (sweep_1) or to the top 


two rows (sweep_2). 


1. General (Point) Manipulation 
The general manipulation of the control points allows the designer to move 
the individual control points in any of the three coordinate directions (r,z,8). Each 


point can be moved in one direction at a time. The axial direction (z) and radial 
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direction (r) are in the same units as are used for the data read by the program. The 
tangential coordinate is in radians. The following Fig.24 shows the effect on the blade 


geometry of moving one (surface edge) control point in the tangential direction. 
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Figure 24. The Effect of Moving Control Point 11 in the Tangential Direction .05 
Radians 


2. Programmed (Shape) Manipulation 

The programmed (shape) manipulations, as currently implemented, move all 
the control points for the six surfaces that share the same row number, simultaneously. 
The tangential movement of the top row or top two rows give different types of blade 
lean. By moving each row individually, different lean schedules can be obtained. 


Figures 25 and 26 show examples of the effect of tangential perturbations. 
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Xe 


Figure 25. Blade Lean Created by Moving the Top Row of Control Points (0.05 
Radians in Each Direction). 


Two special manipulations were implemented for sweep. Both sweep the blade 
tip along the chord direction either forward or backward. The first (sweep_1) moves 
only the top row of control points a prescribed percentage of the rotor chord in the 
required direction (see Fig. 27). The second (sweep_2) moves the top two rows of 
control points simultaneously, the same fraction of chord, in the same way as in 
sweep_l (see Fig. 28). Moving one row of control points will affect the whole blade 
span, but with gradually decreasing effect the further the distance from the row being 
moved. 

The chord slide or ’sweep in the chord direction’ implemented in sweep_l 


and sweep_2 is achieved by a composed (z,#) manipulation. The movement of the 
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Figure 26. Blade Lean Created by Moving the Top Two Rows of Control Points (0.05 


Radians in Each Direction). 


tip section is shown in Fig. 29. The axial displacement required is calculated as 
a percentage of the axial chord. That percentage is the parameter passed to the 
function. The tangential displacement is calculated from the rate of change of the 
tangential coordinate along the chord. Sweep_2 preserves the original blade slope at 


the wall, whereas sweep_l does not. 


D. PROGRAM LANGUAGE AND STRUCTURE 
A listing of the initial version of the code is given in [Ref. 19]. The package 
was implemented using C++ object-oriented programming. The program consists of 


six classes. Some of the classes are used as tools by the main classes. The class is 
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Figure 27. Forward and Aft Sweep (10% of Chord) Using Sweep-l. 


considered an entity with defined functions. The user is allowed only access to class 
members through the class interface (class member functions). This facilitates the 
development because once the class is developed it is not affected by the following 
development in the program. The application that uses the class works through the 
interface and is not affected by the internal implementation. For example, one of 
the member functions in the blade class is the ”read_surf” function that reads the 
pressure and suction surfaces, transforms the data to (z,y,z) coordinates, fits the 
two Bezier surfaces, and then adds the leading and trailing edge surfaces. 

The classes used in the program are named point, vector, matrix, surface, 


bsurface and blade. A brief explanation of each of the classes will follow. 
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Figure 28. Forward and Aft Sweep (10% of Chord Using Sweep_2. 


1. *point’ 

The point class represents a point in space. The class members are three 
variables of type double (z,y,z). Since a point in space is considered a vector from 
the origin to that point, the point class will have some of the vector operations 


implemented. The functions defined for the class are 


1. scalar multiplication and division; 
2. addition, subtraction, inner product and vector norm; 


3. print, single coordinate retrieval (e.g. value of z only). 
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Figure 29. Sweep with Both Forward and Aft-Sweep of 10% of Chord 


2. ‘vector’ 

The vector class represents an array of points of type double. The class was 
constructed to facilitate the matrix vector operations. The class members are the 
vector dimension of type integer and the vector starting address. The functions 


defined for the vector are 


1. vector - constant addition, multiplication, division and subtraction; 
2. vector - vector addition, subtraction, multiplication (cross product); 
3. vector - matrix multiplication; 


4. dimension retrieval, vector norm, vector norm squared, index of maximum 
value, assign value to an element of the vector, element retrieval and range of 
element retrieval. 


dd 


3. ’*matrix’ 
Matrix class members are double array of type double, row dimension of type 


integer and column dimension of type integer. The functions defined for the class are 


1. matrix - matrix multiplication, addition, subtraction; 

2. matrix - vector multiplication; 

3. retrieval of element, submatrix, row dimension, column dimension; 
4. QR factorization; 


5. set the value of an element and resize the matrix dimension. 


4. ’surface’ 
Surface class members are a double array of type point, row dimension and 
column dimension. It is equivalent to a three dimensional array of size row xX col x 3. 


The functions defined for the class are 
1. retrieval of row dimension, column dimension, element and component of an 
element; 
2. surface - matrix multiplication; 
3. surface - surface subtraction; 


4. set value of surface elements, value of individual element with a point, value 
of an individual element with three values of type double. 
5. ‘*bsurface’ 
The bsurface class represents the Bezier control points for a surface. The class 
members are surface of dimension 4 x 4, matrix u for the parameter u, matrix v for 
the parameter v, column dimension and row dimension for the two matrices (both 


have the same dimension). The functions defined for the class are 


1. retrieval of an element, print, surface calculation, one point calculation for 
either an index of the parameters the value of the parameters, resize the pa- 
rameter matrices; 
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2. fit a Bezier surface to given surface, first derivative w.r.t to u or v at each 
point or specific point; 


3. find the parameter u for a given v, z, the parameter v for a given u, z and both 
the parameters u,v for a given 2, z; 


4. calculate conical sections on the surface. 


6. ’blade’ 


The blade class represents the blade geometry with its six bsurfaces, number of 
blade sections, number of points on each of the pressure and suction surfaces, number 


of points on trailing edge and number of points on the leading edge. The functions 


defined for the blade are 


1. read surface data in either data points format or control point format; 


2. print output in PLOT3D format, control points or general format in the (r, z, 0) 
coordinates for grid generation; 


3. control point manipulation and subdivision. 
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Vv. APPLICATION TO THE 
INCORPORATION OF SWEEP 


The geometry package was used to investigate swept geometries. Sweep per- 
turbations were implemented using the package predefined routines. After creating 
a geometry, the blade coordinates were calculated on conical sections for the grid 
generator. ‘The grid generator used was ‘T'CGRID supplied to the Naval Postgradu- 
ate School by the NASA Lewis Research Center [Ref. 20]. The generated grid was 
checked before the CFD solver was applied. The viscous code used here was RVC3D 
’Rotor Viscous Code 3-D”, also supplied by NASA Lewis [Ref. 18]. The code outputs 
the solution in PLOT3D format. However, the solution was examined using FAST, 
"Flow Analysis Software Toolkit”, also from NASA. The geometry was also imported 
into IDEAS through the IGES file written by the package. A solid blade was created 
and the levels of stress were checked for the purpose of comparing the effect of the 
specific geometry manipulations on the stress levels. 

The geometry perturbations included swept forward and swept backward, with 
different values of sweep using sweep_1 and sweep.2. The base geometry used was 


the transonic rotor designed by Sanger [Ref. 7]. 


A. SOFTWARE TOOLS 
1. TCGRID 


The grids were generated using TCGRID (Turbomachinery C-GRID) [Ref. 
20]. The code reads the different formats used by NASA (MERIDL, Crouse-Tweedt, 
or user supplied). The program adds leading- and trailing-edge circles to supplied 
formats except the general format (user supplied). A meridional grid is generated 
between hub and tip, onto which the blade geometry is interpolated. The three- 
dimensional grid is generated from 2-D blade-to-blade grids stacked radially. The 
grid geometry is written to a file in PLOT3D format. 
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2. RVC3D 

RVC3D is a thin-layer Navier-Stokes code with implicit residual smoothing 
[Ref. 20, 18, 21, 22]. The code can be applied as an Euler solver or with either 
Baldwin-Lomax or Cebeci-Smith turbulence models. The equations are formulated 
in cartesian coordinates with rotation then transformed to body-fitted coordinates. 
The viscous terms are neglected streamwise, but retained in the other directions. The 
equations are discretized using a Runge-Kutta scheme. The flow field condition is set 
by specifying the inlet stagnation temperature and pressure, exit hub static pressure 
and the rotational speed. All other parameters will adjust to satisfy the equations 
(inlet axial Mach, exit flow angle etc.). The convergence history 1s printed out and 
the solution is considered converged when the residuals are reduced three orders of 
magnitude and the exit total pressure and total temperature remain nearly constant. 


The solution is written to a PLOT3D-compatible file. 


3. FAST 
FAST was used to examine the solution by displaying different flow field pa- 


rameters calculated from the ”q” vector [Ref. 21]. The program reads in the grid file 
and the solution file associated with it. The program has various capabilities that 
include the generation of contour plots of a function, isosurfaces, grid handling and 


checking, etc. 


4. IDEAS 

IDEAS MASTER SERIES was used to examine the importability of the ge- 
ometry into a CAD/CAM program and to examine the levels of maximum stress 
with the different geometry configurations. The program can read the geometry in 
different formats, manipulate the geometry, and create solids. The stress analysis is 
carried out within the program using the simulation module in which the mesh is 
generated. The blade was clamped at the root section and only the stresses due to 


rotation were examined. 
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B. 


Nelson Sanger [Ref. 7] for the Naval Postgraduate School Turbopropulsion Laboratory 
(TPL). The Sanger stage was a replacement to the transonic compressor designed by 
Vavra for the TPL research test rig [Ref. 23]. 
maximum rpm were dictated by the test rig. The rotor design parameters are shown 


in Table III. 


REFERENCE AND SWEPT GEOMETRIES 


The selected reference geometry was the Sanger rotor which was designed by 


Tip Speed 396.2m/sec 

(1300 ft/sec) 

Design Weight Flow 7.75 kg/sec 
(17.09 lb-m/sec) 


i 
ee 
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Table III. Rotor Design Parameters. 



















The overall rotor dimensions and 


The Sanger rotor was selected as a basis for the implementation of sweep for 


the following reasons: 


. The rotor was designed using CFD. 
. The rotor was aggressively designed, with high loading. 

. The geometry is accessible and the design process is documented. 

. Test results would be available in the near future for validation purposes. 


. If a successful swept design was obtained, it would be possible to build and 


test it in the same rig. 
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The Sanger rotor coordinates were read by the geometry package in MERIDL3 
format and fitted with Bezier surfaces. The control points for the Bezier represen- 
tation are given in Table IV. The leading and trailing edge surface definitions are 
different from the original blade due to the use of the Bezier representation. While 
trying to keep the blade chord as close to the original geometry as possible, with min- 
imum changes in the blade surfaces (i.e. subdivision), the leading surface at the hub 
was found to be blunter than the original blade. While not intended, this did provide 
an opportunity to examine the effect of leading edge variation on the flow field and 
performance. To isolate the effect of sweep, the fitted geometry (’base’ geometry) will 
be used as a reference from which to determine the effect of sweep, rather than the 
original geometry (Sanger geometry). Figures 30 and 31 show the grids for both the 
Sanger and Base geometries. 

The examination of sweep involved several steps, each yielding insight as the 
geometry package was interfaced with the software tools. First, sweep of 5% in 
both forward and backward directions was applied to the geometry using sweep_2. 
The grids for the two swept geometries and the base geometry were generated using 
the same parameters used for generating the Sanger rotor grid. The grid generator 
(TCGRID) was used as supplied without any modification. After examining the 
overall performance computed using RVC3D at near stall conditions, it was seen that 
the base geometry suffered a 0.4% drop in energy-averaged efficiency compared to 
the Sanger geometry. The opposite was expected due to the smoothness conditions 
imposed on the base geometry by the Bezier package. It was found that the handling 
of the leading edge in TCGRID was the reason. An example of the grid at the 
leading edge at mid span is shown in Fig. 32 (The data obtained with this grid 
are referred to as Data Set 1). The use of TCGRID as supplied was clearly not 
suitable because the program recomputes the leading edge points according to a user- 
supplied parameter that controls the spacing of the grid points at the leading edge. 


That spacing is also varied along the blade span linearly relative to the spacing at 
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Table IV. Bezier Representation of Sanger Rotor [Base Geometry]. 
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Figure 30. Grid for Sanger Blade Row 
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Figure 31. Grid for Fitted (Base) Blade Row 
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the hub by another specified parameter. While changing these parameters could 
have improved the leading edge grid spacing at some locations radially, it was an 
undesirable approach due to the shape change that accompanies the gridding process, 
and the lack of control over the other radial stations. Therefore the program was 
changed to compute internally the grid spacing parameter at each station, while 
neglecting the user-supplied parameters. 

The grid then obtained at the mid-span leading edge is shown in Fig. 33, and 
the results obtained with the grid are refered to as Data Set 2. Four cases of sweep 
were examined in addition to the base geometry and the Sanger geometry, all at near 
design conditions. The four cases were swept with 10% tip chord movement along the 
chord direction using sweep_1 and sweep_2 in both forward and backward directions. 
Figures 34, 35, 36 and 37 show the grids generated for the swept blade rows. After 
computing the flow field, the base geometry was still showing lower performance. It 
was not clear at first that the effect was again due to the grid, since it appeared that 
the leading edge bluntness was causing an over-expansion and creating a stronger 
shock. Being unable to explain the increase in relative Mach number with reduced 
mass flow rate, the grid at the leading edge was magnified and was found to be 
relatively coarse. This was due to the use of the same number of grid points over 
the leading edge while having a leading edge of longer arc length. The problem was 
overcome by increasing the number of ‘leading edge’ points supplied to TCGRID and 
decreasing the number of surface points by the same amount. That caused TCGRID 
to consider a portion of the leading edge as a part of the blade surface, giving better 
leading-edge resolution. The final grid is illustrated for the mid-span section in Fig. 38 
and the data obtained with this grid are referred to as Data Set 3. 

Data Set 3 is discussed in the most detail but Data Set 2 is also reported 
(Appendix C) because there are lessons to be learned from the results, and because 
they show that some of the trends determined for the effects of sweep were not grid 


dependent. 
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Figure 33. Grid for Date Set 2. 
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Figure 34. Swept Forward +10% Using Sweep-l 
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Figure 35. Swept Back -10% Using Sweep-1 
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Figure 37. Swept Back -10% Using Sweep_2 
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C. RESULTS OF FLOW ANALYSIS 
iP Overall Performance 


The test cases were examined numerically using RVC3D. The overall ’energy- 
averaged’ performance obtained for the different geometries and leading edge grids is 
summarized in Table V. The exit-static to inlet-total pressure ratio and rotational 
speed for Data Sets 2 and 3 were the same as the design point for the Sanger rotor. 
Data Set 1 was obtained with an exit-static pressure to inlet-total pressure ratio that 
is different from cases 2 and 3. All solutions were obtained at 3000 iterations. The 
RMS density residuals were down to the same level. All cases achieved a reduction 
of over three orders of magnitude. The convergence history is shown in Figure 39 for 
Data Set 3. The back-swept rotor was observed to have the highest efficiency in all 
cases, and to pump a higher mass flow rate. The forward-swept rotor had the lowest 
efficiency, and pumped the lowest flow rate. The Sanger geometry performed better 
than the fitted (base) geometry in Data Set 1 and 2, and the effect of backward sweep 
was to recover that performance. For Data Set 3 the opposite was true and the fitted 
geometry performed slightly better than the Sanger geometry. The effect of sweep on 


performance for Data Set 3 may be deduced from an examination of the flow field. 


2. Exit Flow Field (Peripherally-Averaged) 


The radial distribution of circumferentially energy-averaged efficiency is shown 
for Data Set 3 in Fig. 40. The forward-swept geometry is seen to have lower efficiency 
over the blade span than the base geometry, the Sanger geometry and the swept-back 
geometry. The swept-back geometry showed higher efficiency compared to the base 
geometry. The swept-back rotor efficiency was high from hub to 50% span, and near 
the tip. The efficiency for the base geometry was insignificantly lower than the Sanger 
geometry over a large portion of the blade span, but slightly higher over the outer 
20%. 

Figures 41 - 44 show the radial distributions of exit total pressure to reference 


pressure, exit static pressure to reference pressure, exit total temperature to reference 


of 
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Figure 39. Convergence History for Data Set 3. 
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Case Mass Flow | Efficiency | Total Pressure | Total Temperature 
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Table V. Overall Performance Results. 


temperature and exit static temperature to reference temperature, respectively. The 
total pressure for the swept-forward rotor at the hub is seen to be higher than the 
swept-back rotor. That difference is reduced gradually to zero at about 25% of span, 
staying small to nearly 80% of span, and then reversing in direction. The static 
pressure is higher for the swept-forward and lower for the swept-back compared to 
the base. The total temperature for the swept-forward rotor is higher at the hub 
than for the swept-back rotor. At the tip there is no significant difference between 
the energy added by the different rotors. 

Figures 45 - 18 show exit Mach number, meridional flow angle ¢(arctan(*)), 
total velocity to reference (sonic speed) and swirl angle a(arctan(*#)) where u is the 
axial velocity. The velocity was nearly identical for the four rotors. There is seen to 


be a slight change in the exit flow angles between the rotors near the tip. 
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Figure 40. Radial Distribution of the Circumferentially Energy-Averaged Efficiency. 
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Figure 41. Radial Distribution of the Circumferentially Energy-Averaged Total-to- 
Reference Pressure. 
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Figure 42. Radial Distribution of the Circumferentially Energy-Averaged Static-to- 
Reference Pressure. 
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Figure 43. Radial Distribution of the Circumferentially Energy-Averaged Total-to- 
Reference Temperature. 
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Figure 44. Radial Distribution of the Circumferentially Energy-Averaged Static-to- 
Reference Temperature. 
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Figure 45. Radial Distribution of the Mach Number. 


62 


0.46 


0.44 
Sanger 
: Base 
0.42 Swept-—Forward 
‘1 Swept-Back 


0.4 


Radius (ft) 
Oo 
w 
co 


Oo 
Ww 
oO 


0.34 


0.32 





-100 -80 ~-60 —40 -20 0 20 40 60 80 100 
Phi 


Figure 46. Radial Distribution of the Flow Angle (¢) in the Meridional Plane. 
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Figure 47. Radial Distribution of the Total Velocity to Reference Sonic Velocity. 
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Figure 48. Radial Distribution of the Swirl Angle (a). 


64 


3. Blade Passage Flow Field 

In the following discussion, only the results from Data Set 3 are included. 
Results from Data Set 2 are given in Appendix C however, and based on the overall 
performance and the peripherally-averaged rotor exit profiles, there was consistency 
observed in the effect of sweeping the base geometry either forward or backward. 

The flow field computed for the Sanger rotor at design point (Data Set 3) is 
shown in Fig. 49 - 54. The Sanger rotor Mach number contour plots show that there 
is a mild shock wave at about M = 1.2 at the tip section. The lower half of the 
blade from hub to 50% span has subsonic flow, except for the small region around the 
leading edge. The boundary layer over the suction surface remains relatively thin in 
the flow direction over most of the blade span. Over the outer 10% of the span the 
boundary layer thickens progressively toward the tip. That boundary layer growth is 
affected by the occurrence of the passage shock which started at 75% of span. Also, 
the blade wake at the trailing edge is thicker toward the tip. 

The flow field for the base geometry 1s compared with the flow field for the 
Sanger rotor in Fig. 55 - 60 (Data Set 3). The plots are composite, showing the Mach 
number variation across the blade channels for the two blades side-by-side, with the 
base geometry blade above the Sanger blade. The flow field for the base geometry was 
slightly affected by the leading edge shape. A slight over-expansion over the leading 
edge occurred due to the leading edge geometry. The boundary layer thickness was 


relatively thinner than the Sanger rotor and the wake width was smaller. 
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Figure 50. Mach Number in the Blade Passage of the Sanger Rotor at 25% Span. 
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Figure 52. Mach Number in the Blade Passage of the Sanger Rotor at 75% Span. 
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Figure 54. Mach Number in the Blade Passage of the Sanger Rotor in the Case Wall 
Boundary Layer. 
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Figure 55. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 10% Span. 
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Figure 56. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 25% Span. 





1.494 
: 1.448 
0.803 


0.457 





0.111 


Figure 57. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 50% Span. 
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Figure 58. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 75% Span. 
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Figure 59. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 90% Span. 
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Figure 60. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) in the Case Wall Boundary Layer. 
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Similar composite contour plots of Mach number are shown for the base geome- 
try and swept-forward rotor in Fig. 61 - 66, and for the base geometry and swept-back 
rotor in Fig. 67 - 72. In general, the shock location appears to move forward with 
respect to the leading edge in the swept-forward rotor, and backward with respect to 
the leading edge in the swept-back rotor. 

Figure 73 shows the suction surface pressure and Fig. 78 shows the velocity 
vectors at the third grid point away from the suction surface for the swept-back blade. 
Corresponding plots are given in Figs. 75 and 74 for the base geometry, and in Figs. 
77 and 76 for the swept-forward blade. Examination shows that the location of the 
shock with respect to the blade leading edge at the tip section is moved very little. 
However, the structure of the shock over the blade surface is clearly changed. The 
forward sweep results in a single, continuous shock wave. The backward sweep results 
in two separate shocks. 

The secondary flow created by the shock structure can be seen in the velocity 
vector plots. The radial motion is seen to be largest in the swept-forward rotor and 


least in the swept-back rotor. 
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Figure 61. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Forward Rotor (Upper Blade) at 10% Span. 
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Figure 62. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Forward Rotor (Upper Blade) at 25% Span. 
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Figure 63. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Forward Rotor (Upper Blade) at 50% Span. 
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Figure 64. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Forward Rotor (Upper Blade) at 75% Span. 
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Figure 65. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Forward Rotor (Upper Blade) at 90% Span. 
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Figure 66. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Forward Rotor (Upper Blade) in the Case Wall Boundary Layer. 
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Figure 67. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Back Rotor (Upper Blade) at 10% Span. 
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Figure 68. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Back Rotor (Upper Blade) at 25% Span. 
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Figure 69. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Back Rotor (Upper Blade) at 50% Span. 
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Figure 70. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Back Rotor (Upper Blade) at 75% Span. 
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Figure 71. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Back Rotor (Upper Blade) at 90% Span. 
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Figure 72. Mach Number in the Blade Passage of the Base Rotor (Lower Blade) and 
Swept-Back Rotor (Upper Blade) in the Case Wall Boundary Layer. 
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Figure 73. Pressure on the Suction Surface of the Swept-Back Rotor Blade. 





Figure 74. Velocity Vector Near Suction Surface (3"¢ Grid Point) of the Swept-Back 
Rotor Blade. 
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Figure 75. Pressure on the Suction Surface of the Base Rotor Blade. 
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Figure 76. Velocity Vector Near Suction Surface (3’? Grid Point) of the Base Rotor 
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Figure 77. Pressure on the Suction Surface of the Swept-Forward Rotor Blade. 
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Figure 78. Velocity Vector Near Suction Surface (3°? Grid Point) of the Swept- 


Forward Blade. 
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D. RESULTS OF STRESS ANALYSIS 


The stress analysis for the three geometries in Data Set 3 were carried out 
using IDEAS. The geometry surfaces were transformed to a solid for the the creation 
of the finite element model. The stresses were calculated due to centrifugal loads at a 
rotational speed of 30,000 rpm. The design speed is 27,000 rpm. The swept-forward 
rotor showed an increase in the stress levels on the pressure side to 47,400 psi from 
a stress level of 34,400 psi for the base rotor. The swept-backward rotor showed a 
stress level of 41,700 psi on the suction side. The stress level should be increased for 
the base rotor and swept-forward rotor when the aerodynamic loads are added. For 
the swept-back rotor the stress levels will be reduced due to the aerodynamic load. 


The following figures show the calculated stress levels for the three rotors. 
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Figure 79. Stress Levels for the Base Rotor. 
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Figure 80. Stress Levels for the Swept-Backward Rotor. 
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Figure 81. Stress Levels for the Swept-Forward Rotor. 
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VI. DISCUSSION 


The 3D Bezier-surface geometry representation was successfully implemented 
to represent and manipulate an axial compressor blade geometry. The representation, 
since it depends on fewer parameters than traditional blade-element representations, 
limits the design space. However, it limits the space to acceptable geometries by not 
allowing manipulations that would result in radial waviness and to shapes that are 
inherently machinable. The fewer number of points required for blade representation 
simplifies the process of interfacing aerodynamic and structural analyses, and offers 
promise of optimization. 

The application of the geometry representation to an existing geometry (the 
Sanger rotor) resulted in an improvement in performance. Although the improvement 
was small it was viewed in a positive way since one would expect a poorer performance 
with a blunter leading edge. The improvement might be due to several reasons, all 
of which are equally possible. It could be the minor change in shape introduced in 
the fitting process, or it could be the improved blade surface quality enforced by the 
package. 

The swept cases revealed several important results related to sweep. It was 
shown that the swept-back rotor efficiency and mass flow rate were higher than the 
swept-forward rotor. The shock structure was seen to be weaker for the swept-back 
rotor while stronger for the swept-forward rotor. The results obtained were in agree- 
ment with the reported results in the open literature except for the effect of sweep on 
radial flow reported by Wadia et al [Ref. 13]. They reported that, in the back-swept 
rotor, there was significant radial flow in the boundary layer which weakens the flow 
at the tip section. In the present study the opposite was found. One thing to bear in 
mind is that the sweep schedules are very different and there could be effects other 
than sweep in Wadia et al’s design to cause the difference. They also reported that 


the back sweep reduced the stall margin significantly and reasoned that this was due 


87 


to the boundary layer being centrifuged toward the tip section. In the work reported 
here an examination of stall margin was not conducted; but if the conclusion drawn 
by Wadia concerning the negative effect of centrifuged boundary layer flow is correct, 
we can cautiously expect an improvement in the stall margin using swept-back design 
because a reduction in that effect with backward sweep was observed here. 

In spite of the trouble experienced with the leading edge grid generation it 
was very useful to show that the leading edge could have a dramatic effect on the 
blade row performance and flow structure. Figures 55 - 60 for Data Set 3 can be 
compared with the corresponding composite figures given in Appendix C for Data 
Set 2 to see the strong effect that the grid generation at the leading edge had on 
the predicted flow field. However, it is also interesting to note that the performance 
trends predicted to result from forward or backward sweep were the same for both 
data sets, and therefore are grid independent. 

The difficulty with the grid generation arose because the Bezier representation 
was new, and the leading edge was a larger fraction of the blade chord. It is noted that 
the grid generation that was used was one that was developed for the Q3D design and 
analysis approach. It is possible that a grid generation process could be developed 
that would make use of the analytic properties of the Bezier representation to more 
optimally distribute and orient grid points. 

A final comment concerns the enormous simplicity provided by the Bezier 
representation when going from aerodynamic analysis to stress analysis. The complete 
blade geometry (for the base rotor) is contained in the data in Table IV. (Compare 
this to the MERIDL format in Appendix A). Within IDEAS, the geometry can be 
manipulated and returned (through the geometry package program) to begin a new 
flow analysis cycle. There is a direct and exact interface between the aerodynamic 
and the structural analysis processes, which again points to the need to improve the 


coupling of the geometry package with grid generation for CFD. 
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WALT: CONCLUSIONS AND 
RECOMMENDATIONS 


CONCLUSIONS 


From the presented study the following was concluded: 


. A 3D axial blade row geometry package was implemented successfully using a 
3D Bezier surface representation. 


. The representation is highly flexible and allows manipulation easily and quickly. 


. When the package is fully-developed, MDDO (Multi-Disciplinary Design Op- 
timization) is thought to be feasible due to the fewer number of control pa- 
rameters, compatability with structural analysis and CAD/CAM programs, 
and the expected increase in computer speeds. 


. Sweep has a significant effect on aerodynamic performance, and the new geom- 
etry package will enable a more complete study to be conducted. The present 
results are not sufficient to draw general conclusions concerning the relative 
benefit of forward and backward sweep. 


. Leading-edge gridding has a significant effect on the computed performance 
and flow structure. 


. The package allows different leading-edge configurations to be examined more 
easily than is possible with current design systems. 


RECOMMENDATIONS 


The following recommendations are made: 


. Calculations are required within the geometry package to derive parameters 
that are traditionally used by designers in current quasi-3D-based design pro- 
cesses. 


. A 3D grid generator which is fully compatible with the Bezier representation 
is required. 


. A technique to start a new design using the package should be developed. 


. A graphical user interface should be developed. 
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5. The process of optimization, involving aerodynamic performance, structural 
weight and manufacturing cost, should be examined. 
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APPENDIX A. SANGER ROTOR GEOMETRY 
IN MERIDL3 FORMAT 


MERIDL3 format 1s read into TCGRID as the number of points defining the 


hub and tip followed by the hub and tip coordinate points. Then the blade surface 


coordinates are read into TCGRID after reading in the number of radial sections 


on the blade, the number of points for each section and the number of blades. The 


sequence in which the points are read is, z coordinate for all data points followed 


by y coordinate then the z coordinate. The following is a listing of the Sanger-rotor 


blade definition in MERIDL3 format supplied to TCGRID. 
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APPENDIX B. GENERATING IGES FILES 


Initial Graphics Exchange Specification (IGES) [Ref. 24] is the standard for 
defining geometric objects to ensure importability and system independence. The 
IGES files that are written according to the standard can be read by CAD programs, 


finite element programs, and CNC machines programs. 


1 IGES FILE STRUCTURE 

IGES files have a fixed 80-column ASCII file format. Columns 1-72 include 
the data while columns 73-80 represent the sequence field that numbers the lines 
contained in each of the different sections. The file consists of five sections to define 
the geometry. Each geometry element is considered as a separate entity that has to 
be assigned a pre-defined entity number. The five sections are start, global, directory 


entry, parameter data and terminate. 


a. Start Section 
Start section may contain text in columns 1-72 and a letter ’S’ in column 73. 


One at least is required even if empty, except for the sequence field. 


b. Global Section 

This section contains information describing the program that generated the 
file, and information needed by the reading program. The records are identified with 
the letter ’G’ in colimn 73 and sequenced in column 74-80. The parameters in the 
global section include system ID, precision magnitude, precision significance, scale, 


units, time stamp and version, among others. 


C; Directory Section 
Each entity in the file should have one directory record in the directory entity 
section. Each record is 20 fields of fixed length. Each field is eight characters and the 


filed is written on two lines. The records are identified by the letter ’D’. 
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d. Parameter Data Section 

This section defines the entities to be drawn according to the definition given 
in the IGES standard. Each geometric object (line, circle, etc.) is assigned an "entity 
number” which is followed by a list of numbers separated by commas. The definition 
for the geometric entity can take as many lines as necessary to have a complete 


representation. The records are identified by the letter ’P’. 


e. Terminate Section 

Terminate section is the last line in the file and is always one line. It has 10 
fields of fixed length of eight characters. The fields contain the total number of lines 
of each section. The total number of lines in each of the start section, global section, 
directory entry section and parameter data section are written in the first, second, 


third and fourth fields respectively. The record is identified by the letter ’T’. 


2. GEOMETRY IN IGES FORMAT 


The natural B-spline surface entity (128) is used to represent the blade sur- 
faces. Since the Bezier surfaces are a special case of the B-spline surfaces, the blade 
surfaces are fully represented by this entity. The full blade is represented by six 
B-Spline surface (128) entities [Ref. 24]. 
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APPENDIX C. RESULTS OF FLOW ANALYSIS 
WITH DIFFERENT GRIDS 


The following results are from Data Set 2. 
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Figure 82. Radial Distribution of the Circumferentially Energy-Averaged Efficiency 
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Figure 83. Radial Distribution of the Circumferentially Energy-Averaged Total-to- 


Reference Pressure. 
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Figure 84. Radial Distribution of the Circumferentially Energy-Averaged Static-to- 
Reference Pressure. 
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Figure 85. Radial Distribution of the Circumferentially Energy-Averaged Total-to- 
Reference Temperature. 
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Figure 86. Radial Distribution of the Circumferentially Energy-Averaged Static-to- 
Reference Temperature. 
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Figure 88. Radial Distribution of the Flow Angle (¢) in the Meridional Plane. 
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Figure 89. Radial Distribution of the Total Velocity to Reference Sonic Velocity. 
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Figure 90. Radial Distribution of the Swirl Angle (a). 
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Figure 91. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor 
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Figure 92. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 50% Span. 
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Figure 93. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) at 75% Span. 
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Figure 94. Mach Number in the Blade Passage of the Sanger Rotor (Lower Blade) 
and Base Rotor (Upper Blade) in the Case Wall Boundary Layer. 
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